***Masuoka, Grose and Junn "Sexual Harassment and Candidate Evaluation" Replication code; August 2021

**Replication file for Study 1 from main text
***Variable Recodes
label define pretrtA 1 "Dem" 2 "Rep" 3 "No pty"
label values Pre_treat pretrtA

label define postrtA 1 "control" 2 "placebo" 3 "treatment"
label values Post_treat postrtA

gen treatgrp=0
replace treatgrp=1 if Pre_treat==1 & Post_treat==1
replace treatgrp=2 if Pre_treat==1 & Post_treat==2
replace treatgrp=3 if Pre_treat==1 & Post_treat==3
replace treatgrp=4 if Pre_treat==2 & Post_treat==1
replace treatgrp=5 if Pre_treat==2 & Post_treat==2
replace treatgrp=6 if Pre_treat==2 & Post_treat==3
replace treatgrp=7 if Pre_treat==3 & Post_treat==1
replace treatgrp=8 if Pre_treat==3 & Post_treat==2
replace treatgrp=9 if Pre_treat==3 & Post_treat==3
label define trtgpA 1 "D_control" 2 "D_Placebo" 3 "D_treat" 4 "R_control" 5 "R_Placebo" 6 "R_treat" 7 "NP_Control" 8 "NP_Placebo" 9 "NP_treat"
label values treatgrp trtgpA

gen treatment=0
replace treatment=1 if Post_treat==3

gen female=Q4
recode female 1=0 2=1
rename Q5 educ
rename Q6 income
rename Q3 age
rename Q8 religatt
gen evangel=Q9
recode evangel 1=1 2=0
gen employed=Q10
recode employed 13=1 14=0
gen voted=Q17
recode voted 1=1 2/3=0
gen authr1=Q60_1
recode authr1 1=0 2=1
gen authr2=Q61_1
recode authr2 1=1 2=0
gen authindex=authr1+authr2
gen gendrole=Q18_1+Q18_2+Q18_3
gen gendrole2=gendrole
recode gendrole2 3=1 4=2 5=3 6=4 7=5 8=6 9=7 10=8 11=9 12=10 13=11 14=12 15=13
gen sexism=Q18_4+Q18_5+Q18_6+Q18_7
gen sexism2=sexism
recode sexism2 4=1 5=2 6=3 7=4 8=5 9=6 10=7 11=8 12=9 13=10 14=11 15=12 16=13 17=14 18=15 19=16 20=17
gen feminist=Q19
recode feminist 1=4 2=3 3=2 4=1 
rename Q20_1 therm_fem
rename Q20_2 therm_Dem
rename Q20_3 therm_Rep
gen dem=0
replace dem=1 if Q12==2
gen repub=0
replace repub=1 if Q12==1
gen indep=0
replace indep=1 if Q12==3
gen ptyid=Q12
recode ptyid 4=.
label define partyA 1 "Republican" 2 "Democrat" 3 "Indep"
label values ptyid partyA
gen pre_vote=99
replace pre_vote=1 if Q21DemFpre==5
replace pre_vote=2 if Q21DemFpre==4
replace pre_vote=3 if Q21DemFpre==3
replace pre_vote=4 if Q21DemFpre==2
replace pre_vote=5 if Q21DemFpre==1
replace pre_vote=1 if Q21RepFpre==5
replace pre_vote=2 if Q21RepFpre==4
replace pre_vote=3 if Q21RepFpre==3
replace pre_vote=4 if Q21RepFpre==2
replace pre_vote=5 if Q21RepFpre==1
replace pre_vote=1 if Q21NopartyFpre==5
replace pre_vote=2 if Q21NopartyFpre==4
replace pre_vote=3 if Q21NopartyFpre==3
replace pre_vote=4 if Q21NopartyFpre==2
replace pre_vote=5 if Q21NopartyFpre==1
rename Q24combined pre_therm
gen post_vote=Q26_combined
recode post_vote 1=5 2=4 3=3 4=2 5=1
rename Q29 post_therm
gen diff_vote=post_vote-pre_vote
gen diff_therm=post_therm-pre_therm
gen gendpty=0
replace gendpty=1 if Q12==1 & female==1
replace gendpty=2 if Q12==1 & female==0
replace gendpty=3 if Q12==2 & female==1
replace gendpty=4 if Q12==2 & female==0
replace gendpty=5 if Q12==3 & female==1
replace gendpty=6 if Q12==3 & female==0
label define gdpt1 1 "Rep_women" 2 "Rep_men" 3 "Dem_women" 4 "Dem_men" 5 "Indep_women" 6 "Indep_men"
label values gendpty gdpt1
recode gendpty 0=.
gen match=999
replace match=1 if Q12==1 & Pre_treat==2
replace match=1 if Q12c==1 & Pre_treat==2
replace match=1 if Q12==2 & Pre_treat==1
replace match=1 if Q12c==2 & Pre_treat==1
recode match 999=0
gen genptytrt=0
replace genptytrt=1 if treatgrp==3
replace genptytrt=2 if treatgrp==6
replace genptytrt=3 if treatgrp==9
label define gttA 1 "D_treat" 2 "R_treat" 3 "NP_treat"
label values genptytrt gttA
recode genptytrt 0=.
gen SHdummy=Post_treat
recode SHdummy 1=0 2=. 3=1
gen DWRM=99
replace DWRM=1 if gendpty==3
replace DWRM=0 if gendpty==2
recode DWRM 99=.
label define dummyB 1 "Dem_women" 0 "Rep_men"
label values DWRM dummyB
gen DWDM=99
replace DWDM=1 if gendpty==3
replace DWDM=0 if gendpty==4
recode DWDM 99=.
label define dummyC 1 "Dem_women" 0 "Dem_men"
label values DWDM dummyC
gen DWRW=99
replace DWRW=1 if gendpty==3
replace DWRW=0 if gendpty==1
recode DWRW 99=.
label define dummyD 1 "Dem_women" 0 "Rep_women"
label values DWRW dummyD
gen RWRM=99
replace RWRM=1 if gendpty==1
replace RWRM=0 if gendpty==2
recode RWRM 99=.
label define dummyE 1 "Rep_women" 0 "Rep_men"
label values RWRM dummyE
gen RWDM=99
replace RWDM=1 if gendpty==1
replace RWDM=0 if gendpty==4
recode RWDM 99=.
label define dummyF 1 "Rep_women" 0 "Dem_men"
label values RWDM dummyF
gen DMRM=99
replace DMRM=1 if gendpty==4
replace DMRM=0 if gendpty==2
recode DMRM 99=.
label define dummyH 1 "Dem_men" 0 "Rep_men"
label values DMRM dummyH
gen dwrm_inter=DWRM*SHdummy
gen dwdm_inter=DWDM*SHdummy
gen dwrw_inter=DWRW*SHdummy
gen rwrm_inter=RWRM*SHdummy
gen rwdm_inter=RWDM*SHdummy
gen dmrm_inter=DMRM*SHdummy
gen SHparty=treatgrp
recode SHparty 1/2=. 3=1 4/5=. 6=0 7/9=.
gen dwrm_inter2=DWRM*SHparty
gen dwdm_inter2=DWDM*SHparty
gen dwrw_inter2=DWRW*SHparty
gen rwrm_inter2=RWRM*SHparty
gen rwdm_inter2=RWDM*SHparty
gen dmrm_inter2=DMRM*SHparty
gen ctlparty=treatgrp
recode ctlparty 1=1 2/3=. 4=0 5/9=.
gen dwrm_inter5=DWRM*ctlparty
gen dwdm_inter5=DWDM*ctlparty
gen dwrw_inter5=DWRW*ctlparty
gen rwrm_inter5=RWRM*ctlparty
gen rwdm_inter5=RWDM*ctlparty
gen dmrm_inter5=DMRM*ctlparty
gen PBCLdummy=Post_treat
recode PBCL 1=0 2=1 3=.
gen PBSHdummy=Post_treat
recode PBSHdummy 1=. 2=1 3=0
gen dwrm_inter3=DWRM*PBCLdummy
gen dwdm_inter3=DWDM*PBCLdummy
gen dwrw_inter3=DWRW*PBCLdummy
gen rwrm_inter3=RWRM*PBCLdummy
gen rwdm_inter3=RWDM*PBCLdummy
gen dmrm_inter3=DMRM*PBCLdummy
gen dwrm_inter4=DWRM*PBSHdummy
gen dwdm_inter4=DWDM*PBSHdummy
gen dwrw_inter4=DWRW*PBSHdummy
gen rwrm_inter4=RWRM*PBSHdummy
gen rwdm_inter4=RWDM*PBSHdummy
gen dmrm_inter4=DMRM*PBSHdummy
gen repwomen=0
replace repwomen=1 if gendpty==1
gen repmen=0 
replace repmen=1 if gendpty==2
gen demwomen=0
replace demwomen=1 if gendpty==3
gen demmen=0
replace demmen=1 if gendpty==4
gen othergp=1
replace othergp=0 if gendpty==1
replace othergp=0 if gendpty==2
replace othergp=0 if gendpty==3
replace othergp=0 if gendpty==4
gen sameparty=999
replace sameparty=1 if dem==1 & treatgrp==3
replace sameparty=1 if repub==1 & treatgrp==6
replace sameparty=0 if dem==1 & treatgrp==1
replace sameparty=0 if repub==1 & treatgrp==4
recode sameparty 999=.
gen dwrm_inter6=DWRM*sameparty
gen dwdm_inter6=DWDM*sameparty
gen dwrw_inter6=DWRW*sameparty
gen rwrm_inter6=RWRM*sameparty
gen rwdm_inter6=RWDM*sameparty
gen dmrm_inter6=DMRM*sameparty
gen oppparty=999
replace oppparty=1 if dem==1 & treatgrp==6
replace oppparty=1 if repub==1 & treatgrp==3
replace oppparty=0 if dem==1 & treatgrp==4
replace oppparty=0 if repub==1 & treatgrp==1
recode oppparty 999=.
gen dwrm_inter7=DWRM*oppparty
gen dwdm_inter7=DWDM*oppparty
gen dwrw_inter7=DWRW*oppparty
gen rwrm_inter7=RWRM*oppparty
gen rwdm_inter7=RWDM*oppparty
gen dmrm_inter7=DMRM*oppparty
gen samep1=99
replace samep1=1 if Q12==2 & Pre_treat==1
replace samep1=1 if Q12==1 & Pre_treat==2
replace samep1=0 if Q12==2 & Pre_treat==2
replace samep1=0 if Q12==1 & Pre_treat==1
recode samep1 99=.
gen sameparty2=999
replace sameparty2=1 if dem==1 & treatgrp==2
replace sameparty2=1 if repub==1 & treatgrp==5
replace sameparty2=0 if dem==1 & treatgrp==3
replace sameparty2=0 if repub==1 & treatgrp==6
recode sameparty2 999=.
gen dwrm_inter8=DWRM*sameparty2
gen dwdm_inter8=DWDM*sameparty2
gen dwrw_inter8=DWRW*sameparty2
gen rwrm_inter8=RWRM*sameparty2
gen rwdm_inter8=RWDM*sameparty2
gen dmrm_inter8=DMRM*sameparty2
gen sameparty3=999
replace sameparty3=1 if dem==1 & treatgrp==2
replace sameparty3=1 if repub==1 & treatgrp==5
replace sameparty3=0 if dem==1 & treatgrp==1
replace sameparty3=0 if repub==1 & treatgrp==3
recode sameparty3 999=.
gen dwrm_inter9=DWRM*sameparty3
gen dwdm_inter9=DWDM*sameparty3
gen dwrw_inter9=DWRW*sameparty3
gen rwrm_inter9=RWRM*sameparty3
gen rwdm_inter9=RWDM*sameparty3
gen dmrm_inter9=DMRM*sameparty3
gen ctldwdm=999
replace ctldwdm=1 if gendpty==3 & treatgrp==1
replace ctldwdm=0 if gendpty==4 & treatgrp==1
recode ctldwdm 999=.
gen ctldwrw=999
replace ctldwrw=1 if gendpty==3 & treatgrp==1
replace ctldwrw=0 if gendpty==1 & treatgrp==4
recode ctldwrw 999=.
gen ctldwrm=999
replace ctldwrm=1 if gendpty==3 & treatgrp==1
replace ctldwrm=0 if gendpty==2 & treatgrp==4
recode ctldwrm 999=.
gen ctldmrw=999
replace ctldmrw=1 if gendpty==4 & treatgrp==1
replace ctldmrw=0 if gendpty==1 & treatgrp==4
recode ctldmrw 999=.
gen ctldmrm=999
replace ctldmrm=1 if gendpty==4 & treatgrp==1
replace ctldmrm=0 if gendpty==2 & treatgrp==4
recode ctldmrm 999=.
gen ctlrwrm=999
replace ctlrwrm=1 if gendpty==1 & treatgrp==4
replace ctlrwrm=0 if gendpty==2 & treatgrp==4
recode ctlrwrm 999=.
gen ctldwdm2=999
replace ctldwdm2=1 if gendpty==3 & treatgrp==4
replace ctldwdm2=0 if gendpty==4 & treatgrp==4
recode ctldwdm2 999=.
gen ctldwrw2=999
replace ctldwrw2=1 if gendpty==3 & treatgrp==4
replace ctldwrw2=0 if gendpty==1 & treatgrp==1
recode ctldwrw2 999=.
gen ctldwrm2=999
replace ctldwrm2=1 if gendpty==3 & treatgrp==4
replace ctldwrm2=0 if gendpty==2 & treatgrp==1
recode ctldwrm2 999=.
gen ctldmrw2=999
replace ctldmrw2=1 if gendpty==4 & treatgrp==4
replace ctldmrw2=0 if gendpty==1 & treatgrp==1
recode ctldmrw2 999=.
gen ctldmrm2=999
replace ctldmrm2=1 if gendpty==4 & treatgrp==4
replace ctldmrm2=0 if gendpty==2 & treatgrp==1
recode ctldmrm2 999=.
gen ctlrwrm2=999
replace ctlrwrm2=1 if gendpty==1 & treatgrp==1
replace ctlrwrm2=0 if gendpty==2 & treatgrp==1
recode ctlrwrm2 999=.

**Difference in difference analysis
regress post_vote SHdummy DWRM dwrm_inter
regress post_vote SHdummy DWDM dwdm_inter
regress post_vote SHdummy DWRW dwrw_inter
regress post_vote SHdummy RWRM rwrm_inter
regress post_vote SHdummy RWDM rwdm_inter
regress post_vote SHdummy DMRM dmrm_inter

regress post_vote SHparty DWRM dwrm_inter2
regress post_vote SHparty DWDM dwdm_inter2
regress post_vote SHparty DWRW dwrw_inter2
regress post_vote SHparty RWRM rwrm_inter2
regress post_vote SHparty RWDM rwdm_inter2
regress post_vote SHparty DMRM dmrm_inter2

regress post_vote PBCLdummy DWRM dwrm_inter3
regress post_vote PBCLdummy DWDM dwdm_inter3
regress post_vote PBCLdummy DWRW dwrw_inter3
regress post_vote PBCLdummy RWRM rwrm_inter3
regress post_vote PBCLdummy RWDM rwdm_inter3
regress post_vote PBCLdummy DMRM dmrm_inter3

regress post_vote PBSHdummy DWRM dwrm_inter4
regress post_vote PBSHdummy DWDM dwdm_inter4
regress post_vote PBSHdummy DWRW dwrw_inter4
regress post_vote PBSHdummy RWRM rwrm_inter4
regress post_vote PBSHdummy RWDM rwdm_inter4
regress post_vote PBSHdummy DMRM dmrm_inter4

regress post_vote ctlparty DWRM dwrm_inter5
regress post_vote ctlparty DWDM dwdm_inter5
regress post_vote ctlparty DWRW dwrw_inter5
regress post_vote ctlparty RWRM rwrm_inter5
regress post_vote ctlparty RWDM rwdm_inter5
regress post_vote ctlparty DMRM dmrm_inter5

regress post_vote sameparty DWRM dwrm_inter6
regress post_vote sameparty DWDM dwdm_inter6
regress post_vote sameparty DWRW dwrw_inter6
regress post_vote sameparty RWRM rwrm_inter6
regress post_vote sameparty RWDM rwdm_inter6
regress post_vote sameparty DMRM dmrm_inter6

regress post_vote oppparty DWRM dwrm_inter7
regress post_vote oppparty DWDM dwdm_inter7
regress post_vote oppparty DWRW dwrw_inter7
regress post_vote oppparty RWRM rwrm_inter7
regress post_vote oppparty RWDM rwdm_inter7
regress post_vote oppparty DMRM dmrm_inter7

regress post_vote sameparty2 DWRM dwrm_inter8
regress post_vote sameparty2 DWDM dwdm_inter8
regress post_vote sameparty2 DWRW dwrw_inter8
regress post_vote sameparty2 RWRM rwrm_inter8
regress post_vote sameparty2 RWDM rwdm_inter8
regress post_vote sameparty2 DMRM dmrm_inter8

regress post_vote sameparty3 DWRM dwrm_inter9
regress post_vote sameparty3 DWDM dwdm_inter9
regress post_vote sameparty3 DWRW dwrw_inter9
regress post_vote sameparty3 RWRM rwrm_inter9
regress post_vote sameparty3 RWDM rwdm_inter9
regress post_vote sameparty3 DMRM dmrm_inter9

**Differences across treatment groups, accounting for party/gender/party*gender of respondent
tabstat post_vote, stats(mean sd n) by (treatgrp)
anova post_vote treatgrp
tukeyhsd treatgrp

tabstat post_therm, stats(mean sd n) by (treatgrp)
anova post_therm treatgrp
tukeyhsd treatgrp

tabstat diff_vote, stats(mean sd n) by (treatgrp)
anova diff_vote treatgrp
tukeyhsd treatgrp

tabstat diff_therm, stats(mean sd n) by (treatgrp)
anova diff_therm treatgrp
tukeyhsd treatgrp

tabstat post_vote if dem==1 & female==1, stats(mean sd n) by (treatgrp)
anova post_vote treatgrp if dem==1 & female==1
tukeyhsd treatgrp if dem==1 & female==1

tabstat post_vote if dem==1 & female==0, stats(mean sd n) by (treatgrp)
anova post_vote treatgrp if dem==1 & female==0
tukeyhsd treatgrp if dem==1 & female==0

tabstat post_vote if repub==1 & female==1, stats(mean sd n) by (treatgrp)
anova post_vote treatgrp if repub==1 & female==1
tukeyhsd treatgrp if repub==1 & female==1

tabstat post_vote if repub==1 & female==0, stats(mean sd n) by (treatgrp)
anova post_vote treatgrp if repub==1 & female==0
tukeyhsd treatgrp if repub==1 & female==0

tabstat post_vote if Post_treat==3, stats(mean sd n) by (gendpty)
anova post_vote gendpty if Post_treat==3
tukeyhsd gendpty if Post_treat==3

tabstat post_vote if Post_treat==2, stats(mean sd n) by (gendpty)
anova post_vote gendpty if Post_treat==2
tukeyhsd gendpty if Post_treat==2

tabstat post_therm if Post_treat==3, stats(mean sd n) by (gendpty)
anova post_therm gendpty if Post_treat==3
tukeyhsd gendpty if Post_treat==3

tabstat post_therm if Post_treat==2, stats(mean sd n) by (gendpty)
anova post_therm gendpty if Post_treat==2
tukeyhsd gendpty if Post_treat==2

tabstat post_vote if treatgrp==1, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==1
tukeyhsd gendpty if treatgrp==1

tabstat post_vote if treatgrp==2, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==2
tukeyhsd gendpty if treatgrp==2

tabstat post_vote if treatgrp==3, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==3
tukeyhsd gendpty if treatgrp==3

tabstat post_vote if treatgrp==4, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==4
tukeyhsd gendpty if treatgrp==4

tabstat post_vote if treatgrp==5, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==5
tukeyhsd gendpty if treatgrp==5

tabstat post_vote if treatgrp==6, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==6
tukeyhsd gendpty if treatgrp==6

tabstat post_vote if treatgrp==7, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==7
tukeyhsd gendpty if treatgrp==7

tabstat post_vote if treatgrp==8, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==8
tukeyhsd gendpty if treatgrp==8

tabstat post_vote if treatgrp==9, stats(mean sd n) by (gendpty)
anova post_vote gendpty if treatgrp==9
tukeyhsd gendpty if treatgrp==9

tabstat diff_vote if dem==1 & female==1, stats(mean sd n) by (treatgrp)
anova diff_vote treatgrp if dem==1 & female==1
tukeyhsd treatgrp if dem==1 & female==1

tabstat diff_vote if dem==1 & female==0, stats(mean sd n) by (treatgrp)
anova diff_vote treatgrp if dem==1 & female==0
tukeyhsd treatgrp if dem==1 & female==0

tabstat diff_vote if repub==1 & female==1, stats(mean sd n) by (treatgrp)
anova diff_vote treatgrp if repub==1 & female==1
tukeyhsd treatgrp if repub==1 & female==1

tabstat diff_vote if repub==1 & female==0, stats(mean sd n) by (treatgrp)
anova diff_vote treatgrp if repub==1 & female==0
tukeyhsd treatgrp if repub==1 & female==0

tabstat diff_therm if dem==1 & female==1, stats(mean sd n) by (treatgrp)
anova diff_therm treatgrp if dem==1 & female==1
tukeyhsd treatgrp if dem==1 & female==1

tabstat diff_therm if dem==1 & female==0, stats(mean sd n) by (treatgrp)
anova diff_therm treatgrp if dem==1 & female==0
tukeyhsd treatgrp if dem==1 & female==0

tabstat diff_therm if repub==1 & female==1, stats(mean sd n) by (treatgrp)
anova diff_therm treatgrp if repub==1 & female==1
tukeyhsd treatgrp if repub==1 & female==1

tabstat diff_therm if repub==1 & female==0, stats(mean sd n) by (treatgrp)
anova diff_therm treatgrp if repub==1 & female==0
tukeyhsd treatgrp if repub==1 & female==0

******Effect Size analysis
by Post_treat, sort: esize twosample post_vote, by(DWRM) all
esize twosample post_vote if Post_treat==3, by(DWRM) all
esize twosample post_vote if Post_treat==3, by(DWRW) all
esize twosample post_vote if Post_treat==3, by(DWDM) all
esize twosample post_vote if Post_treat==3, by(RWRM) all
esize twosample post_vote if Post_treat==3, by(RWDM) all
esize twosample post_vote if Post_treat==3, by(DMRM) all
by treatgrp, sort: esize twosample post_vote, by(DWRM) all
by treatgrp, sort: esize twosample post_vote, by(DWRW) all
by treatgrp, sort: esize twosample post_vote, by(DWDM) all
by treatgrp, sort: esize twosample post_vote, by(RWRM) all
by treatgrp, sort: esize twosample post_vote, by(RWDM) all
by treatgrp, sort: esize twosample post_vote, by(DMRM) all

**Placebo analysis
tabstat post_vote if gendpty==3, stats(mean sd n) by (Post_treat)
anova post_vote Post_treat if gendpty==3
tukeyhsd Post_treat if gendpty==3

tabstat post_vote if gendpty==1, stats(mean sd n) by (Post_treat)
anova post_vote Post_treat if gendpty==1
tukeyhsd Post_treat if gendpty==1

tabstat post_vote if gendpty==4, stats(mean sd n) by (Post_treat)
anova post_vote Post_treat if gendpty==4
tukeyhsd Post_treat if gendpty==4

tabstat post_vote if gendpty==2, stats(mean sd n) by (Post_treat)
anova post_vote Post_treat if gendpty==2
tukeyhsd Post_treat if gendpty==2

**Feeling thermometer analysis
tabstat post_therm if dem==1 & female==1, stats(mean sd n) by (treatgrp)
tabstat post_therm if dem==1 & female==0, stats(mean sd n) by (treatgrp)
tabstat post_therm if repub==1 & female==1, stats(mean sd n) by (treatgrp)
tabstat post_therm if repub==1 & female==0, stats(mean sd n) by (treatgrp)

**same party analysis (same party candidate, harassment vs control)
regress post_therm sameparty DWRM dwrm_inter6
regress post_therm sameparty DWDM dwdm_inter6
regress post_therm sameparty DWRW dwrw_inter6
regress post_therm sameparty RWRM rwrm_inter6
regress post_therm sameparty RWDM rwdm_inter6
regress post_therm sameparty DMRM dmrm_inter6

**opposite party analysis (harassment vs control)
regress post_therm oppparty DWRM dwrm_inter7
regress post_therm oppparty DWDM dwdm_inter7
regress post_therm oppparty DWRW dwrw_inter7
regress post_therm oppparty RWRM rwrm_inter7
regress post_therm oppparty RWDM rwdm_inter7
regress post_therm oppparty DMRM dmrm_inter7
